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Abstract 

The correlation functions of two or more Euclidean Wilson loops 
of various shapes in Euclidean anti-de Sitter space are computed by 
considering the minimal area surfaces connecting the loops. The sur¬ 
faces are parametrized by Riemann theta functions associated with 
genus three hyperelliptic Riemann surfaces. In the case of two loops, 
the distance L by which they are separated can be adjusted by con¬ 
tinuously varying a specific branch point of the auxiliary Riemann 
surface. When L is much larger than the characteristic size of the 
loops, then the loops are approximately regarded as local operators 
and their correlator as the correlator of two local operators. Similarly, 
when a loop is very small compared to the size of another loop, the 
small loop is considered as a local operator corresponding to a light 
supergravity mode. 
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1 Introduction 


Since the inception of the AdS/CFT correspondence [T], many aspects of 
gauge theory have been made clearer by studying relevant quantities in string 
theory. One important example is the case of the Wilson loop, which is one 
of the most fundamental quantities of a gauge theory. According to HE], a 
Wilson loop can be computed in string theory by considering a minimal area 
surface in anti-de Sitter space that ends on the loop in the boundary of the 
space. 

Many examples of Wilson loops have been studied in this context |1] , with 
much more general examples recently computed in mum where a new, 
inhnite parameter family of minimal area surfaces were studied. In those 
works the minimal area surfaces were constructed analytically in terms of 
Riemann theta functions associated to hyperelliptic Riemann surfaces closely 
following previous work by Pi- 

Our objective in this paper is to holographically compute the correlation 
function of multiple Wilson loops of various shapes by means of minimal 
area surfaces. The AdS/CFT correspondence implies that the minimal area 
surfaces have these loops as boundary. For two correlated Wilson loops, when 
one Wilson loop is “very small” we will associate it with a particle emitted 
by the worldsheet ending on the other loop n. Then the correlator in this 
case may by viewed as the correlator of a Wilson loop and a chiral primary 
operator. Most of the technique used in this paper was established in mu 
and the reader is encouraged to review them. 

In a conformal held theory, it is often very useful to know the correlation 
functions of nonlocal operators like the Wilson loop. In the limit when the 
t’Hooft coupling A is small the Wilson loop may be perturbatively expanded. 
For large A, however, one may compute OPE of a Wilson loop correlator 
either by taking the correlator of the loop with the relevant chiral primary 
operator or by taking the correlator of the loop with another loop |10j . 

For example, the correlation of two circular Wilson loops of equal size 
a separated by a distance L has been well studied where, interestingly, it 
is shown that a phase transition occurs when L is large enough [mini [131 
HD. For small L the surface connecting the two loops is the catenoid. The 
transition occurs when L gets large enough to exceed a critical value when 
the minimal area surface prefers to separate into two hemispheres connected 
by a vanishingly narrow tube. 

The connected part of the correlators of two Wilson loops in the super- 
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gravity limit is 


{W\C^)W{C2))conn = exp(-;^) (1.1) 

/TTQ! 

where S is the Nambu-Goto action 



Additionally, others [T^IT^ITT] have also studied the correlator of null Wilson 
loops. 

The paper is organized as follows: In section two we review how to go 
from the sigma model to a nonlinear scalar equation of motion. Section three 
contains the main technical difference with the previous solutions [U]. We set 
up the auxiliary Riemann surface in section four. Then we use the Riemann 
surface in section hve to compute examples with two loops. In the following 
section we show how it may be extended to more than two Wilson loops, 
then we make concluding remarks in the hnal section. 

2 Sigma Model in Euclidean AdS^ 

A convenient way to imagine Euclidean AdS^ is to consider it as a subspace 
of dehned by the equation 

-Xl + Xl + Xl + Xl = -l, (2.3) 

with an obvious S'0(3,1) = S'L(2,C) global invariance. The action of the 
string in conformal coordinates is given by 

S = {dX^X'^ - A{X^X^^ -1)) da dr 

= ^j^{daXd^X + daYd^Y + daZd^Z)dadT, (2.4) 

where A is a Lagrange multiplier and the /r indices are raised and lowered 
with the TZ^’^ metric 

3 

{v,w) =— VqW^ . (2.5) 

i=l 
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The complex coordinates z and z are related to the world sheet coordinates, 
(Jq = {a,T) hy z = a + ir, and z = a — ir. X,Y, Z are called Poincare 


coordinates (see (2.24)). The string equations of motion are given by 


ddX^ = AX^ , 

where A, the Lagrange multiplier is given by 

A = ax^ax‘‘ = {x,,x,). 


( 2 . 6 ) 


(2.7) 


The conformal condition is encoded in the Virasoro constraint 

(X„X,) = 0 = (X„X,-). (2.8) 

The first step to solving the equations of motion is to reduce the problem 
to an equation with a single unknown scalar field. In Euclidean AdS^ space 
this scalar equation is the cosh-gordon equation 


dda = 4 cosh a . 


(2.9) 


This reduction mechanism referred to as Pohlmeyer reduction [IH] has been 
used in the study of many related problems. In the context of Minkowski 
space-time this procedure was used by Jevicki and Jin [T9| and by Kruczenski 
[20] to hnd new spiky string solutions, by Irrgang and Kruczenski [21] to 
study new open string solutions in AdS^, and by Alday and Maldacena [22| 
to compute certain light-like Wilson loops. In a more geometric guise it was 
employed [8] to study constant mean curvature surfaces in hyperbolic space. 
We review the idea and follow closely what was done there. 

In Euclidean AdS^, form a basis 


where 


f = (W,W„W„7V) 

(X, N) = (X„ N) = (X„ N) = 0, (X, N) = 1 


Note that (2.3) and (2.8) imply 

{X,X,) = {X,X,) = 0 

and 

(X„X,,) = (X„X,,) = (X„X,,) = (X„X,,) = 0 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 
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respectively. Define 


(X,, A'_.) := 2e“, (A„-, N) := 2ff*e“, (A„, iV) := A* 


(2.14) 


where is the Hopf differential and is the mean curvature of the surface 
described by the solution to the string equation of motion. Since we are 
concerned here with a minimal area surface (described by the string equations 
of motion) we will have a vanishing mean curvature and the second equation 
in (2.14) is equal to zero. We want to study what happens when the basis 


undergoes a small motion with the hope that second derivative quantities will 
tell us something useful about the sigma model. For this we write second 
derivatives as a linear combination of the basic vectors. Thus we write 


Xzz = cbX + bXz + cXz + dN. 
Then taking inner product with X gives 

-a = {X,Xz,) = -{Xz,Xz) 


(2.15) 


which implies a = 2e“ due to (2.14). Doing same for the other vectors in the 
basis gives 


0 = {Xz, Xzz) = c{Xz, Xz) ^ c = 0 
0 = (X„ X,,) = 6(X„ X,) ^6 = 0 

0 = {N,Xzz) = -d 


This shows that 

X,, = 2e“X, (2.16) 

which is exactly the equation of motion (2.6) taking A = 2e“. Note that 
without the minimal area condition, i.e. = 0, equation (2.16) will have 
an additional term proportional to H^N. Repeating this for the other second 
derivatives we get 


X,, = «,X, - A^N 


and 


X-z-z = a-zX-z-A^N. 

Consider motion of the basic vectors given by 
ti,z ti,z 


Vijtj . 


(2.17) 
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Note that since second derivatives are expressed entirely in terms of first 
derivatives, quantities such as ti^z and ti^zz are all expressible in terms of 
hrst derivatives. Consequently, the matrices U and V contain a and its 
derivatives, and (and its complex conjugate). Taking second derivatives 
and imposing the compatibility condition ti^zz = ti,zz leads to 

{Uij^z — Vij^z)tj + {UijVjk — VijUjk)tk = 0 . (2.18) 

This equation may be written in matrix form, after dropping the t, as 


Uz-Vz + [U,V] = 0 


(2.19) 


with U and V determined from (2.17) as 


/ 

0 

1 

0 

0 

\ 


/ 

0 

0 

1 0 


0 

az 

0 

A’^ 


, ^ = 


2e“ 

0 

0 0 


2e“ 

0 

0 

0 



0 

0 

Of A’^ 

V 

0 

0 

-\A^e-^ 

0 



V 

0 

-\A^e-^ 0 0 


The compatibility equation implies that is a holomorphic function. Fur¬ 
thermore, it leads to a generalized cosh-gordon equation 


1 


- 2e“ - ^A^A^e-^ = 0 . 


( 2 . 20 ) 


If in (2.20) we scale the Hopf differential A^ —)■ 2^4^ and similarly for its 


conjugate, the generalized cosh-gordon equation becomes 


- 2e“ - 2A^A^e-^ = 0 . 


( 2 . 21 ) 


Changing coordinates to w = z,w = z followed by a transfor¬ 
mation of the held a —>■ d = ck — ^ logA^A^, the resulting equation is the 
standard cosh-gordon equation (2.9). This procedure is equivalent to setting 
A^ = 2 in ( |2.20[ ). The solution for a in terms of theta functions is found to 
be 

a(z,z) = -, (2.22) 


e 


Ai/2 

A2/2 


iC{z,z)) 


where C is a complex vector of dimension equal to g, the genus of the Riemann 
surface, and Ai, A 2 are the characteristics of the theta function. They are 
also vectors of integers with dimension g. 
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To connect the scalar field a with the solutions for (2.6) we introduce a 
hermitian matrix 

■ X0 + X3 X,-tX 2 
Xr+iX 2 X0-X3 
Then construct Poincare coordinates 

1 


X = 


(2.23) 


Z = 


X 


X 


22 


■V ^21 


(2.24) 


The final connection is to recognize that (2.20) is also the compatibility 
condition for the system of equations [S] 


= f/$ 


= 1/$ 


where U and V are given by 




0 2Ae“/2 

2g-«/2 Q, 




2 a/2 

X 


0 


(2.25) 


(2.26) 


The parameter A is the spectral parameter which emerges in the study of 
integrable differential equations. Given a solution for the pair of equations 


(2.25), the Poincare coordinates are related to ( |2.25| ) by 
$11$12 + 4*22^21 


X + iY = 




12^12 


$ 22^22 


Z = 


V det $ det $1 




12^12 


4*22*^ 


(2.27) 


22^22 


where $ = 


One advantage of Poincare coordinates is that 


$11 $12 

^21 ^22 

the boundary of AdS^ is now a copy of TZ^. 

Riemann theta function solution for the cosh-gordon equation was found 


in m, then the system (|2.25[) was solved in terms of those theta functions. 


The solution for $ furnishes solutions for (2.27) for various values of |A| = 1 


In those works the spectral parameter did not coincide with branch points 
of the auxiliary Riemann surface by which the family of Wilson loops are 
parametrized. The Wilson loops are in the boundary of Euclidean AdS^ 
where Z = 0. 


3 New Solutions 

Technically, the major deviation from the previous work [6] is that we set one 
branch point equal to A. We need |A| = 1 and throughout we will use A = 1. 
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In order to explore the behavior of the solutions under the restriction that A 
and a branch point are equal, we begin by reviewing the previous solutions. 
The solutions for (2.27) was given in equations (3.94) and (3.95) of [5] as 


Z = 


X + iY = e 


+0 

moiomi 

^flZ+UZ 

2 


|Ac-/”)I^ + |9(C-/”)!' 


(3.28) 


2fs+2i>. - Ip‘WC + J") - ^(C - Jf'WC + Ip ‘) p 


i9(c-/;;)p + ie(c-/r)i 


with 


CPi 


fi = -Dp^\ne{ ), 


u = —Dp^ In 6 


Pi 


'pi 


Ai/2 

A2/2 


rP4 


(3.30) 


'pi 


Pi and Pi are points on the Riemann surface. Concretely, pi = 0 and p 4 = A, 
where the value of A is taken to lie on the unit circle. 

When, as in the present case, A coincides with a branch point, Z vanishes 
and X + iY may be singular. The automatic vanishing of Z can occur if 
6{2 /g~^), where x is a branch point equal to one, is automatically zero. 
This problem can be traced back to the fact that at A = x the columns of <I> 


fail to be unequal and the system (2.25) becomes singular in the sense that 
det < 1 ) = 0 . 

In order to fix this we must regularize $ and choose a new matrix con¬ 
sisting of a column from the old < 1 > and a column from the regularized ^reg- 


The solution for the system (2.25) is 


$ = 


/ e(c+A{\))e(d) 0{c-A{x))0id) -i;{\,z,z) \ 

0iO0{d+A{X)) 0{O0id-A{X)) 


^(C+d(A)) 0{d) i/)(A,z,z) _ ^(C~d(A)) 0{d) —'ip{X,z,z) 


.(3.31) 


\ 0{C)0{d+A{X)) 


0{C)0id-AiX)) 


/ 


n is a function on the Riemann surface equal to y/X. i4(A) = is the Abel 
map and d = i{c, c, c^)'^ with c G C"’, c° G M is an arbitrary column vector. 6 
is a shorthand notation for a theta function with characteristics Ai, A 2 , and 
= z dQoo + ^ /oo ^^0 where (Aloo,o are the Abelian differentials 


of the second kind . They simplify to (3.30) when x 7 ^ A = 1. 





















We introduce a local parameter q = y/\ — x in the vicinity of y and take 


derivatives dq\q=o of (2.25) . Since dq commutes with dz,z we get 

dzdq^q=, = dqU\q=,^ + Udq^q=o 
dzdq^q=0 = dqV\q=0^ + Vdq^q=0 


(3.32) 


The matrices U and V depend on A only through A itself, and since dqA = 
dq{q^ + x)) therefore, dqlJ and dqV both vanish at g = 0. So the upshot is 
that if <h satisfy (2.25) then dq^\q=Q also does. 

The elements of <hreg := dq^ are[^ 


$ 


reg,ll 


^reg,21 
^reg,12 
^reg,22 


0(0 

1 

W) 

V 

W) 


diO 


{dq[9{C + A(A))et’(^’^>")] - e{C + Al(A)) dq ln0(d + Al(A))} 

{^^[^(C + Al(A))e^(^’^’"i] - e{C + Al(A)) dq ln0(d + Al(A))} 

{dq[9{C - Al(A))e-^(^’"’")] - 9{C - yl(A)) dq \n9{d - Al(A))} 


{aj0(C - Al(A))e-^(^’"’")] - 9{C - Al(A)) dq \n9{d - Al(A))} 


(3.33) 


where A = + y. 


Two solutions for (2.25) with non zero determinants can now be con¬ 


structed from <h and ^reg by taking the following combinations; 


$1 = 


^ 0 (d-A(A))^ 

^reg,ll --*^*12 


<h 

'if r 


m 

e{d) ®22 


*ho — 


e(d+A(A)) ^ $ 


eid) 


reg,12 


e{d+A{X)) ^ $ 


. (3.34) 


■reg,21 —'=t'22 J \ “'±'21 'i^reg,22 

In this paper we focus on the solution <I)i. The new solutions in terms of 


Poincare coordinates can be derived using (2.27) as 


X + iY = {a;i(A, z, z) [dq'ip{X, z, z) — dq ln6*(y4(A) -|- d)] + X 2 (A, z, z)} 

(3.35) 


Z = 


\v\-^ e{c-A{\)) e(c-A(A))+e(c-A(A)) ejc-Ajx)) 

Her 


^An overall factor of which does not depend on C has been dropped. This 

amounts to dropping a constant multiplier in the determinant of $reg and therefore does 
not affect the solutions in any significant way. 
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where the functions Xi and X 2 are 


xi(A,C) = 


|n|2 0(C - A(A)) e{C + A(A)) - 0(C - Al(A)) 0(C + A(A)) 


0(C - ^(A)) e{C - Al(A)) + 0(C - Al(A)) 0(C - A{\)) 


X2(A,C) = 


|n|^ 0(C - AI(A)) d,e{c + AI(A)) - 0(C - AI(A)) d,e{c + ai(a)) 


0 (C - Al(A)) 0(C - ^(A)) + e{C - A(A)) 0(C - A{X)) 


The determinant det $1 is given by 

det $1 = K 


^^(C) 


(3.36) 


(3.37) 


It may be shown that the constant k is equal to 

—V 


K = 


e 


[dgO 


Ai/ 2 
Ao/2 
Ai/2 
Ao/2 


( 0)0 


Ao/2 

A 2/2 


(A(A)0 


( 0 ) 

Ao/2 ■ 
A 2/2 


(A(A)) + 9 


Ai/2 

Ao/2 


{A{X))d,9 


Ao/2 

A2/2 


ag0(d - A{x))e{d + A{x)) + e{d - A{x))dge{d + a{x))- 
e{d + A{x))e{d - A{x)) - 


One important thing to note here is that upon replacing det $1 in (3.35) by 
(3.37), it becomes clear that the zeros of Z correspond to the zeros of 9{Q 
and 0(C)- In this paper we focus on the Wilson loops that are generated by 
the zeros of 0(C)- 

The function C’ is bilinear in ^ and A and is given by; 


'0(A, z, z) = Di In 


i( 0 )_ 

«(h) 


i-z + z). 


(3.39) 


where Di In 


m 

SiloA 


is a real constant. Thus, we have that dq'ip, which appears 


in ^reg, IS given by 


dq'lp = dqdijj {-z + z) . 


(3.40) 


Note that . With these expressions, namely (3.39) and (3.40), 

it can be checked that <I>i ($ 2 ) satishes the pair of equations in (2.25) and 
the solutions (|3.35|) indeed describe the strings moving in AdS^. 


(A(A)) 


(3.38) 
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Figure 1: g=3 hyperelliptic Riemann surface along with a choice of basis 
cycles. 

4 The Auxiliary Riemann Surface 

It is important to setup the Riemann surfaces which parametrize our solu¬ 
tions. In particular we will show that moving one particular branch point 
around changes the separation between two correlated Wilson loops. Please 
refer to m [911231 EH Ea Eg EZ] for a review of hyperelliptic Riemann surfaces 
and how they parametrize our solutions. 

Consider a genus g = 3 Riemann surface which is given by a hyperelliptic 
curve dehned by 

= A(A-a)(A + l/a)(A-6)(A-6)(A+c)(A + c), aeR,b^ceC (4.41) 

The corresponding hyperelliptic Riemann surface with a canonical basis of 
cycles is displayed in Figure^ Due to the symmetry A —?• — 1/A, the branch 
point c is related to 6 by c = —1/6. This leaves two free parameters, a 
and b, that determine the hyperelliptic Riemann surface. Of these two we 
hx a = X = 1 and study when A = x- So given b we can compute (, and 
the Riemann period matrix which leads to a solution for a. This solution 
for a then leads to a solution for the string equations of motion. We can 
change the values of the imaginary or the real part of b and observe the 
behavior of the Wilson loops that emerge. The real part of b is restricted by 
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Refc] 



— Re[b] = 0.9 

— Re[b] = 0.8 
Re[b] = 0.7 
Re[b] = 0.6 

— Re[b] = 0.5 

— Re[b] = 0.4 

— Re[b] = 0.3 

— Re[b] = 0.2 


Figure 2: Several curves for each value of Re{b) arranged in descending order. 


0 < Re{b) < 1 and the imaginary part is constrained by the condition that 
Re{c) = Re{—l/b) < —1. These restrictions can be deduced by inspection 
of Figure^^ 

So given a value for 0 < Re{b) < 1 we can select an appropriate Im{b) 
and study the solutions we obtain. Since Re{c) is a smooth function of Im{b), 
we can plot several curves each for a different value of Re{b) as in Figure^ 
Each curve describes a set of Wilson loops and as we move along the curve 
we observe a change in the space of correlated Wilson loops. 

In order to illustrate this, we study in detail some Wilson loops corre¬ 
sponding to Re{b) = 0.9, 0.6 and 0.4. 


5 Correlation of Two Wilson Loops 


Let us consider the simpler case of the two Wilson loops separated by a 
distance L. We set Re{b) = 0.9 and observe that as we dial the Im{b) from, 
say, 0.002 to 0.009, the separation L between the two Wilson loops decreases 
until they start to intersect. In order words one Wilson loop is moving toward 
the center of the other, and their correlation is expected to become stronger. 
This may be conhrmed by comparing the computed regularized areas for 


different values of Im(b) and taking into account that (1.1) has a negative 
sign in the exponential function (see ra6/e[^. 

If the characteristic size of the Wilson loops is and their separation is 
L, we are going to consider hrst the case L ^ Ob and show how the distance 
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L can be changed by changing the positions of the branch point. Afterwards 
we consider the case where one loop shrinks to a point. 


5.1 Explicit Results for Big Loops Correlator with In¬ 
termediate L 

Let us hrst work out an example when L is of the order of the size of the 
Wilson loops. The results are explicit so that the procedure may be easily 
repeated by any interested person and for other values of the parameter b. 
We can determine the characteristics Ai, A 2 from 


1 /.OO 1 1 

—n = / oj = “A2 T —B.^i 

2 Jo 2 ^ 2 ' 

and inspection of Figure as 



■ 1 ■ 


1 

Ai = 

0 

, A 2 — 

1 


1 

0 


1 


(5.42) 


(5.43) 


The constant vector d may be taken for example to be d = i(l + i, 1 ~ L 1)^ 
but for simplicity we use d = ( 0 , 0 , 0 )^ throughout our calculations. 

The parameter b is chosen from the top most curve in Figure corre¬ 
sponding to Re{b) = 0.9. We start with a small value Im{b) = 0.002 for the 
imaginary part of 6 , and compute (^ = U z + iV z to be; 


C{z,z) = 


■ 0.3814 ■ 


■ 0.3814 ■ 

0.0372 

;2 + 

0.1091 

0.1091 


0.0372 


z. 


(5.44) 


The period matrices denoted by A and B are easily computed to be 


2.4868 i -6.1192 -7.5547 
A= I 1.2817i -6.7992 6.7992 
-2.4868 i -7.5547 -6.1192 


(5.45) 


B = 


0.8250 

-0.5 + 1.8750 i 
0.0141 i 


i 

0.8250 

-0.8250 
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-0.8250 
0.0141 i 
0.5 + 1.8750 i 


(5.46) 










Note that as a Riemann period matrix, R is a 3 x 3 symmetric matrix with 
a positive dehnite imaginary part. 

The matrix 

1 0 0 \ 

0 0 1 (5.47) 

0 10 / 

allows us to compute the complex conjugate of theta functions as 



e 


a 


{C) = e 


—T ■ a 
T-/3 


(T-O 


with T ■ B ■ T = —B. Notice that = 1 and that 


(5.48) 


iT-U = V, 


(5.49) 


which implies that T ■ ( = Furthermore, we have that T ■ Ai 2 = — Ai 2 


The purpose of this excursion is to say that (5.47), (5.48), and (5.49), along 


with the evenness of 6{() and the oddness of 6{() implies that both 6{() and 
9{C) are real. Finally it can be checked that is identically zero. 

Let’s introduce the coordinates of the string worldsheet z = a + ir, 
z = a — IT. We plot the curves describing the zero locus of 9{Q as in Figure 
. Since the Wilson loops are described by the zeros of Z in the boundary 
of EAdSs, the pair of open curves therefore correspond to the Wilson loops 
which are the boundary of the minimal area surfaces in anti-de Sitter space. 
The region lying between the pair of curves corresponds to the minimal area 
surface. Using +iY{Q and Z(() the curves and the region 

lying between them are mapped to Euclidean anti - de Sitter space as shown 
in Figure From the hgure it is clear that one feature of the minimal area 
surface we hnd here is that it self intersects. This is made more explicit in 
where lines of constant sigma are depicted. This feature is present 
in all the examples we will compute in the rest of this paper. 

Another way to think about these solutions is to take the solution for 
the concentric Wilson loops in Euclidean anti - de Sitter space [6] and slowly 
move the inner loop to a region outside the bigger loop. The donut-like 
surface would deform into something similar to the results here. However 
as we will show subsequently, many other higher order correlators can be 
computed here that are not related to our previous results in this simple 
way. 
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Figure 3: Two blue curves describing the zeros of 9{() in worldsheet coor¬ 
dinates for b = 0.9 + 0.002i. The Wilson loops on the right correspond to 
the pair of blue curves and the region between the curves is the minimal area 
surface that connects the two loops. 
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Figure 4: Minimal area surface conn 


ecting the loops in Figure]^ 


Y 



Figure 5: Self intersecting surface connects the loops. 
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5.2 Big Loops Correlator with Small L 


Now that the procedure for our study of Wilson loop correlator has been 
made explicitly clear in the example for intermediate L, we will present the 
results for the example when b = 0.9 + 0.007h Remember that the loops 
approach each other as we increase the imaginary part of b from 0.002 to 
0.009. For values well beyond that point we begin to see the emergence of 
multiple loops from which one can compute the correlator of more than two 
Wilson loops. This pattern of transition from correlator of two loops to more 
than two loops is present in all the lines corresponding to the various real 
part of b plotted in Figure However, these transitions occur along the 
various curves for different values of the imaginary part of b. 

One motivation for computing the case of Im{b) = 0.007 is to show how 
the results for the same pair of open curves in Figure compares when the 
Im{b) is increased. We know that as the loops get closer this indicates a 
stronger correlation, so we expect to see an increase in the absolute value of 
the finite part of the area. 

Again the hyperelliptic Riemann surface is as in Figure and this leaves 
us with the same theta characteristic as before. In this case the corresponding 
pair of open curves and Wilson loops are shown in Figure and the minimal 
surface in Figure\^ 

We have a full analytic result, (9.68) and (9.76), for computing the reg¬ 
ularized area for the minimal area surface connecting the Wilson loops. We 
applied Stoke’s Theorem in separating the hnite part of the area from the 
divergent part and showed that the divergent part is the length of the Wilson 
loop (see the Appendix). Numerically the total area is computed at several 
different values of e using the formula A = 4 f da dre'^ where a = 21og|. 
We then fit the data to the linear model eA = e A finite + L and find both 
the regularized area and the length of the Wilson loops. In Table we show 
these results, confirming that the closer the Wilson loops the more strongly 
correlated they are. 


5.3 Small Loops Correlator 

We now present the results for the case of small Wilson loops in which one 
loop shrinks to a point. An example of such a case is shown in Figure\^foT 
b = 0.4 -|- 0.002 7 If one loop is very small compared to the other, the small 
loop may be viewed as a local chiral operator. This operator corresponds to 
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Figure 6: Two blue curves describing the zeros of 6{() in worldsheet coor¬ 
dinates for b = 0.9 + 0.007i. The Wilson loops on the right correspond to 
the pair of blue curves and the region between the curves is the minimal area 
surface that connects the two loops. 
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0,5 z 
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Figure 7: Minimal area surface connecting the loops in Figure p 
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Table 1: Area of connecting surfaces computed using both numerical and 
analytical methods 


b 

Area by numerical method 

Area by (9.68) and (9.76) 

0.9 + 0.007i 

16.3934 + 

e 

16.3956 + 

€ 

0.9 + 0.002i 

16.3852 + ^9-3767 

e 

16.3853 + ^^-3767 
€ 


X 



Figure 8: Minimal surface connecting a small loop with another loop that has 
shrunk to a point, b = 0.4 + 0.002 T 

a particle emitted by the worldsheet that ends on the big loop |T0] . Another 
possibility is when both loops are small and the distance separating them 
L is much larger than their characteristic size, a^. In this situation one can 
imagine the loops being approximated by local operators with their correlator 
corresponding to the correlator of two-point functions. This interpretation 
can also be applied to cases with more than two Wilson loops. The area of 
the surface connecting the small loop to the point is 

6.1328 

A = -5.2012 + —-— . (5.50) 


19 

















6 Correlator of Three or More Loops 

Now that we have a clear understanding of how the correlation function for 
two Wilson loops in the supergravity limit is computed, the idea may be 
easily extended to higher number of correlators. Recall that in the two loops 
case the Wilson loops and their dual minimal surfaces arose from a pair of 
open curves bounding a horizontal strip in the worldsheet coordinates. It is 
natural, therefore, to expect that higher number of loops will be obtained 
from inserting closed curves between a pair of open curves. So for a correlator 
of three Wilson loops we insert one closed curve, for a correlator of four 
Wilson loops we insert two closed curves, and so on. The more insertions 
we add the more complicated the surface becomes. In some instances some 
of the corresponding Wilson loops may intersect, or may become concentric. 
In our previous work [6] all the Wilson loops obtained were either concentric 
or, in the case of more than two loops, self-intersecting. The present work is 
an extension of those results to nonintersecting loops. We study an example 
with three nonintersecting loops. 

Consider the case when b = 0.6 -|- 0.15 z which correspond to the Re{b) = 
0.6 curve in Figure\^ This setup gives correlators of three Wilson loops since 
the open curves have a single closed curve in between see Figure In this 
example we have one large Wilson loop and two relatively small loops. Using 
the analytic formulae, the area of the surface connecting the three loops is 

26.7055 

71 =-10.1663 +-^-. (6.51) 

It might be interesting to consider a limit where the two small loops 
become point-like. 

7 Conclusion 

We have holographically computed examples of correlators of two or more Eu¬ 
clidean Wilson loops. In the case of two loops we show that their separation 
may be adjusted by changing the location of the branch cut of the auxiliary 
Riemann surface. In some examples where one or more loops may shrink to a 
point, such small loops indicate local chiral operators corresponding to light 
supergravity modes. In previous work [6] we provided solutions of minimal 
surfaces ending on several intersecting loops and in view of that the current 
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Figure 9: Two open curves and a closed curve describing the zeros of 
in worldsheet coordinates for b = 0.6 + 0.15i. The small loops on the right 
correspond to the pair of open curves and the large loop to the closed curve. 
The region between the pair of open curves and the closed curve corresponds 
to the minimal area surface that connects the three loops. 
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Figure 10: Minimal area surface connecting the three loops in Figure 9 


results may be regarded as an extension of those results to nonconcentric and 
nonintersecting loops of various shapes. 
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9 Appendix 

The analytic formulas used in this paper to calculate the area of the minimal 
surfaces are derived in this appendix. 
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9.0.1 Area 


The necessary formulae for computing the regularized area for the minimal 
area surfaces have already been developed and used in [5l [6]. We review 
them here briefly. In the current context, the area allows us to calculate 
the correlation of several Wilson loops or the amplitude of exchange of light 
supergravity modes emitted by the worldsheet ending on a loop. 

Recall that the area of the minimal area surface is given by 

A = 4 e^'dadr, (9.52) 


which diverges at the boundary due to the vanishing of 0{(), and therefore 
requires a regularization According to the AdS/CFT prescription one 
should cut the surface a.t Z = e and write the area as 

A = ^ + Af, (9.53) 

where L is the total length of the Wilson loops and Aj is the finite part which 
is identified with the correlator of the Wilson loops through: 

(Wt(C'i)W(C'2))eonn = , (9.54) 


where A is the ’t Hooft coupling of the gauge theory (not to be confused 
with the spectral parameter). This prescription is equivalent to subtracting 
the area A = ^ of a string ending on the contour of length L and stretching 
along Z from the boundary to the horizon. 

Using Fay’s trisecant identity, we find an expression for the exponential in 


(9.52) to be a sum of a finite term and a term that diverges at the boundary 


where 9 = 0; 


e“ = In 0(0) - In 0(C) 

= iAp3p3 ln0(O) — 9c)ln0(C). 


(9.55) 

(9.56) 


Integrating the second term at the boundary of the surface obviously 
leads to divergence so we need to regulate it. In order to do that we observe 
we may write Z as a product of a non vanishing function and 0 

Z=l9(()lh(z,z) . (9.57) 

^The integral diverges because a = 2 In |, and 9 is zero on the boundary of the surface. 
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Substituting (9.56) into (9.52) and applying Stoke’s theorem we get the ex¬ 
pression 


A = In 0(0) J dadr + j) h ■'Vlnh di — (p n ■ V \n Z di. (9.58) 
In [5] we showed that this can be expressed as 


A = dDpjpg ln0(O) / dadr + / dadrv'^ In h — / In Z . (9.59) 


The preimage of the Wilson loops are a pair of open curves [possibly 
with closed curves in between] in the world sheet coordinates as shown in 
Figure Therefore, Stoke’s Theorem is an effective way to compute the 
area of the surface bounded by the Wilson loops. This surface corresponds 
to the horizontal strip between the open curves in the worldsheet coordi¬ 
nates. Hence, integrating along the open curves and along the two vertical 
boundaries, one at a = do and the other at a = u/. Here do and cx/ are two 
points along the open curves such that cx/ — cxq is the period of X{a, r), i.e. 
W(ao,r) = X(a/,r). 

We parametrize the sine-like curves by the variable ex, so that the lower 
curve is now given as ri((x) and the upper one by r 2 (a). According to Stoke’s 
Theorem, for any smooth real-valued functions Q and P on a regular domain 
D in we have that 


^ dxdy = [ Pdx + Qdy . 
n^dx dyJ ^ JoD 


(9.60) 


We apply (9.60) the hrst term on the right hand side of (9.59) which we 


denote by Aconst- Taking Q = a/2 and P = —t 12 we get 


A 


const 


= -2D 


PlP3 


In0(0)^ / crT 2 {a)—af / dr— / aT[{a)da— / T 2 {a)da+ / Ti(a)da 


(9.61) 

where the subscripts 1,2,3,4 on the integrals indicate left, top, right, and 
bottom boundaries of the domain in the a — t plane. Of course the simplest 
thing to do here is by following elementary calculus and directly write 


Aconst = 4:Dp^p^ In0(0) / (r 2 (cx) - Ti{a))da 


(9.62) 
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However, we chose to arrive at (9.61) by means of (9.60) because it will turn 


out that this approach is very useful for the other two more complicated 


terms in (9.59) where the condition that leads to (9.62) is absent. 


For the second term in (9.59), denoted by A^, we may take Q = dfjlnh 


and P = — In h. This gives 


Ah — 


dr dfj In h\^ 


da d„\xih\j.^T 2 {<j) + / dr drj\xih\ 




+ / da d„\\ih\^^T[{a) + / dadr\^h\^^ — / dadT-\^h\^^. (9.63) 


Similarly, the last term in (9.59) indicated by A^ becomes. 


A, = 


+ 


dTd„\nZ\^^ - 


dr drj In Z 




— / da da-\yi Z\T- 2 T 2 {a) 


I da dfjln Z\T-iT[{a) + J dadr\nZ\ 


da d-r In Z I 


(9.64) 


Since Z vanishes along the open curves parametrized as Ti{a) and T 2 {a), 
we see clearly that only the terms in parenthesis in ( 9.64[ ) are hnite leav¬ 
ing all integrals along sides 2 and 4 which are the two horizontal open curves 
bounding the world sheet to diverge. This is one nice thing about the Stoke’s 
Theorem approach because the divergent part is exposed in very clear man¬ 
ner. To remedy the divergence, we cut the surface at a height e very close 
to the original boundary, and the integrals are no longer divergent up to the 
boundary of this cut surface. This is why the string theory is said to have 
an infrared divergence, but the corresponding gauge theory has an ultravi¬ 
olet divergence. The preimage of the boundary of the cut surface is then 
parametrized by two horizontal curves, ti(cr), t 2 (cr) that lie very close to the 
original ones and on the inside of the worldsheet. Once, this is done the 
formula then becomes 


41. = - 


dr da InZl 




* 2 ( 0 ) rt 2 {o-f) 

drdalnZlao - / 

tl(0) Jti{af) 

- dadaZlt^t'^ia) - j da daZ\t^t[{a) - j da drZ\t^ + j dadrZ\t^'^ 

(9.65) 

The hrst term in parenthesis above consists of integrals along the left and 
right vertical boundaries of the strip and they vanish. So Az = —A^iv where 
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Adiv is the other grouped item along with its coefficient 1/e. Finally, it is 
clear that the total area of the minimal area surface can be written as the 
sum of a term that is hnite and a term that diverges as 1 /e 


with 

and 


^4 -Aconv -^div 

(9.66) 

A —A _L 4, 

^conv -^const ~ 5 

(9.67) 


Adiv = --([ dadaZlt^t^ia)- [ da d^Z\t^t[{a)- j dadrZ\t^+f dadrZ\t^ 

^ ^ J T2 d T\ d r2 d Tl 

(9.68) 

According to the regularization prescription, the divergent part of the 
total area of the minimal surface should be equal to L/e where L is the 
length of the Wilson loop. This implies that if (9.68) is correct, it should 
give the total length of the Wilson loops. 


Recall from |5] that the length of the Wilson loop is 


L = 


' 2+4 


\vz\dl-- 


' 2+4 


|VZ| 


dl. 


(9.69) 


The regularized area which is the hnite part of the total area of the surface 
ending on the Wilson loops is obtained by 

A finite ~ ^ ~ ~ ■ (9.70) 

On the other hand we have 

Az= log ZdadT= / V log Z-dl (9-71) 

J D JdD 


where dl = fidl with n the outward normal vector. With the tangent vector 
to the curve given by {da, dr) we take dl = {dr,—da). Going around the 
loop as before, we obtain 


71 . 




(9.72) 
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We have seen that the hrst term in parenthesis vanishes. Cutting the surface 
at Z = e, leaves us with the relation 


~Az 


\yZ\dl, at Z = e . 


' 2+4 


Hence, the hrst term in (9.69) is exactly equal to the eA^iv found in (9.68). 
So the formula for the regularized area of the surface ending on the Wilson 
loop becomes 


A finite A^^nst "h A^ -|- 


' 2+4 




dl 


(9.73) 


Or more explicitly, 


A finite = 4£>pip3 log^(O) / “ ri{a))da+ j Vlogh ■ d/+- 




dl 


IdD 


12+4 |VZ| 

. (9.74) 

From Z = 6 h we can compute that at Z = 0 we h ave VZ = V6 h and 
V^Z = v'^6 h + 2 ■ Vh. When substituted into (9.74) we get 

A finite = 4F>pip3 log6'(0) / (r2(cr) -Ti{a))da+ [ ^ ■ dl 


>dD 


1 

2 


' 2+4 




dl -\- 


^dl (9.75) 


'2+4 |V0| h 


Looking at the formula for VZ = V6' h it is clear that VZ and V6 are in the 
same direction so that the unit normal may be taken to be —V6/\V6\. This 
further simplihes the above equation for A finite giving an expression purely 
in terms of theta functions and the parametric curves Ti and + 2 ; 


A 


finite 


= 4T)p,p3 log 0(0) j (+2(a) - Ti{a))da + 
= —29'|i7i3 log0(O) y zdz + ^ 




dl 


't2{(j) Jrx(a) 


T2((t) 

Dx \ogO{Qdz 

(9.76) 


The equivalence of lines one and two of (9.76) is shown in [6]. 
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